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ABSTRACT 

We present the first characterization of the excess continuum emission of accreting T Tauri stars 
between optical and near-infrared wavelengths. With nearly simultaneous spectra from 0.48 to 2.4 fiin 
acquired with HIRES and NIRSPEC on Keck and SpeX on the IRTF, we find significant excess 
continuum emission throughout this region, including the /, Y, and J bands, which are usually 
thought to diagnose primarily photospheric emission. The lYJ excess correlates with the excess in 
the V band, attributed to accretion shocks in the photosphere, and the excess in the K band, attributed 
to dust in the inner disk near the dust sublimation radius, but it is too large to be an extension of the 
excess from these sources. The spectrum of the excess emission is broad and featureless, suggestive 
of blackbody radiation with a temperature between 2200 and 5000 K. The luminosity of the lYJ 
excess is comparable to the accretion luminosity inferred from modeling the blue and ultraviolet 
excess emission and may require reassessment of disk accretion rates. The source of the lYJ excess is 
unclear. In stars of low accretion rate, the size of the emitting region is consistent with cooler material 
surrounding small hot accretion spots in the photosphere. However, for stars with high accretion rates, 
the projected area is comparable to or exceeds that of the stellar surface. We suggest that at least 
some of the lYJ excess emission arises in the dust-free gas inside the dust sublimation radius in the 
disk. 

Subject headings: accretion, accretion disks — protoplanetary disks — stars: formation — stars: 
pre-main-sequence 



1. INTRODUCTION 

Classical T Tauri stars (CTTS) are low-mass stars in 
the final stages of disk accretion. Accretion from the 
disk to the star is thought to be governed by the stel- 
lar magnetosphere, where a sufficiently strong magnetic 
field truncates the disk at several stellar radii and guides 
infalling material along field lines to the stellar surface 
at high latitudes, terminating in accretion shocks. His- 
torically, a major factor contributing to the understand- 

j ing that T Tauri stars are accreting matter from their 
circumstellar disks was the spectrum of their excess con- 
tinuum emissi on from ultraviolet to radio wavelengths 

I (jBertoutlll989( ). While the disks themselves are the pri- 
mary source of continuum radiation at wavelengths long- 
ward of 2 fim, optical and ultraviolet continuum emission 
in excess of the photosphere is attributed to accretion 
shocks on the stellar surface. Spectral lines also played a 
key role in developing the current magnetospheric accre- 
tion paradigm, providing the kinematic evidence that the 
stellar magnetosphere channels accreting material from 
the disk to the star and that accretion-powered w inds 
arise from near-stellar regions (|Naiita et al.|[2000l ). In 
this paper we focus on the excess continuum emission in 
a region that has not received close scrutiny to date - 
the spectral region between the optical and the near in- 
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frared, where the photospheric emission from these late- 
type stars peaks. 

The well studied optical and ultraviolet excesses are 
used to derive a fundamental parameter of T Tauri stars 
- the rate of accretion from the disk to the star. The 
inferred accretion rates span several orders of magni- 
tude, from extremes of 10^^° to 10~^ Mq yr~^, with 
a median of ~ 10~^ M(7) yr^-^ at 1 Myr t hat declines 
with increasing age (jHartmann et al.l 119981) . The most 
reliable estimates for disk accretion rates come from 
comparisons of spectrophotometric observations over a 
broad wavelength range to shock models. This was done 
successfully by Calvet and Gullbring (1998, hereafter 
0098;), who accounted for the excess continuum between 
0.32 and 0.52 /im with shock models that attribute the 
Paschen continuuum to optically thick post-shock gas 
in the heated photosphere and the Balmer continuum 
and Balmer jump to optically thin gas in the pre-shock 
and attenuated post-shock regions. This combination of 
spectrophotometry with shock models gives a fairly uni- 
form result for the spectrum of the excess in the Paschen 
continuum, with temperatures ~ 6000 — 8000 K and a 
shape that is essentially blackbody. A more commonly 
used approach to determine accretion rates is to eval- 
uate the excess emission over a limited range of wave- 
length in the Paschen continuum by comparing the depth 
of photospheric absorption lines to those of a template 
matched in temperature, gravity, and projected rota- 
tional velocity. In the presence of a continuum excess, 
photospheric features will be weakened, and a quantity 
known as veiling^ defined as the flux ratio r = ^excess / F^, , 
can be derived (jBasri fc Batalhal 119901 : iHartigan et al.l 
|1991|) . Accretion rates are then determined from a bolo- 
metric correction based on an isothermal slab model 
at an assumed temperature, density, and optical depth 
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(jValenti et al.l ll993t iHartigan et aLlll995t GuUbri ng et 
al. 1998a, hereafter! GHBC98I : IHartigan fc Kenvod [20031: 
IHerczeg fc Hillenbrandll2008l) . 

At near-infrared wavelengths, CTTS show excess emis- 
sion from 2 to 5 fim that is well described by ~ 
1400 K blackbody radiation and is attributed to a raised 
rim of dust at the dust sublimation radius in the in- 
ner disk (Muzero l le et al. 2003, hereafter IMCHD03I : 
iFolha fc EmCTsm] [200H IJohns-Kruh fc Valentil l2Q0l . 
The magnitude of the excess is proportional to the ac- 
cretion luminosity, requiring that the dust be heated by 
radiation from both the photosphere and the accretion 
shock. K-hand interferometry, which loca tes dust in a 
ring a t a few tenths of an AU fro m the star (|Eisner et al.l 
I2005t iMillan-Gabet et al.l [20071) . provides further evi- 
dence that 2 to 5 fim emission in CTTS comes from 
sublimating dust in the inner disk. 

Evidence has been accumulating that there is an unex- 
plained source of excess emission in CTTS between the 
optical and near infrared that cannot be attributed to 
shocks at the base of magnetic funnel flo ws or dust sub- 
limating in the disk. As early as 1990, iBasri fc Batalhal 
found that the veiling longward of 0.5 /xm does not 
steadily decline with wavelength as expected from hot 
accretion shock models but instead is relatively constant 
between 0.5 and 0.8 ^m, and high ve iling at 0.80-0.85 iim 
has m ore recently been reported bv Hart igan fc KenvonI 
(120031) and White & Hillenbrand (2004). Furthermore, 
[Edwards^ al.i (2006) found high veiling at 1 /im, which 
is unlikely to come from 1400 K dust in the inner disk 
since this emission will fall rapidly to short wavelengths 
from its peak around 3 ^m. Also, K-hand interferometric 
studies of T Tauri stars find the angular size of the near- 
infrared emission in modeling of multiple-baseline obser- 
vations (Akeson et al. 2005) and the size- w avelength be- 
havio r in spectrally dispersed observations ([Eisner et al.l 
l2009f) to suggest that gaseous material inside the dust 
sublimation radius is present at a temperature higher 
than that of the dust. However, the lack of a systematic 
study of the wavelength dependence of excess emission 
between 0.5 and 2 /im has inhibited the understanding 
of its role in accreting systems. 

The presence of unexplained continuum emission in 
CTTS between 0.5 and 2 fj,m not only offers the oppor- 
tunity to improve our understanding of the structure of 
accretion disk systems; it also poses a practical problem. 
If we do not understand the shape of the veiling spectrum 
at wavelengths longer than 0.5 /im, then we cannot reli- 
ably convert veiling measurements at these wavelengths 
via a simple bolometric correction into accurate accretion 
luminosities and disk accretion rates. There may also be 
ramifications for extinction determinations of CTTS if 
they are based on far-red colors assumed to be predomi- 
nantly photospheric. 

In this paper we determine the excess spectrum over 
a broad wavelength region between 0.48 and 2.4 /im for 
a sample of 16 classical T Tauri stars in Taurus. The 
presentation includes § 2 describing the sample and data 
reduction, § 3 describing the derivation of the veiling 
and the excess emission spectra, § 4 characterizing the 
behavior of the excess emission spectra, and § 5 describ- 
ing simple models for the excess, followed by a discussion 
in § 6 and conclusions in § 7. We demonstrate that the 
two traditional sources of optical and near-infrared ex- 



cess emission, a hot component arising from an accretion 
shock-heated photosphere with small filling factor and 
a cool component arising from the dust sublimation ra- 
dius of the accretion disk, are not sufficient to describe 
the observed excess. A third component of intermediate 
temperature seems to be required to fit the excess emis- 
sion in this region, with a luminosity of the same order 
of magnitude as derived from the hot shock-heated gas. 

2. SAMPLE AND DATA REDUCTION 

Our sample consists of the 16 classical T Tauri stars 
in Table [T] They were chosen to have a broad range of 
published mass accretion rates, covering four orders of 
magn itude. All were observed with SpeX ([Ravner et al.l 
I2003D at the Infrared Telescope Facility (IRTF) on 2006 
November 26 and 27. A few days later, on November 30 
and December 1, most were also observed nearly simul- 
taneously with HIRES (Vogt et al. 1994) on Keck I and 
NIRSPEC ([McLean et anil998l ) on Keck II. We combine 
the three data sets to address the excess emission over a 
broad wavelength range. 

Seven of the objects in the sample are members of 
binary systems. All but one of the sy stems have sep- 
arations > 0.88" ([White fc Ghezll200l( ). and they were 
resolved by all three instruments. In these cases we ob- 
served only the primary, rotating the spectrograph slit if 
necessary to avoid contamination by the secondary. They 
are noted with an "A" in Table [TJ DF Tan A and B, on 
the other hand, were unresolved at a separation of only 
0.09", and thus both components contribute to the spec- 
tra. The components have similar s pectral types (M2.0 
and M2. 5: 'Hartiga n fc Kenvon| | 2003f) and a flux ratio at 
K of 1.62 (White fc Ghedl2001h . In this work we treat 
DF Tau as a single star with the published properties of 
the primary. 

2.1. SpeX 

The SpeX data were taken in the short-wavelength 
cross-dispersed (SXD) mode, featuring a 0.3" x 15" slit 
and yielding spectra that extend from 0.8 to 2.4 iim at a 
resolving power R = 2000. The detector is a 1024 x 1024- 
pixel Aladdin InSb array with a 0.15" pixel scale. Total 
exposure times ranged from 16 to 48 minutes for the pro- 
gram stars, with 8.2 < J < 10.7, yielding a continuum 
S/N > 250 at J, increasing to longer wavelengths. 

To facilitate the subtraction of sky emission lines, the 
observing sequence for a given object consisted of multi- 
ple 120 s exposures. After the flrst and third exposures 
of a four-exposure sequence, the telescope was nodded 
by 8" (an ABBA pattern) such that the object remained 
in the slit. Normal AO stars near the targets and at sim- 
ilar airmass (usually a difference of less than 0.05) were 
observed to allow the removal of atmospheric absorption 
lines and provide an approximate flux cal ibration. 

Th e data were reduced with Spextool ([Gushing et al.l 
|2004|) . an IDL package containing routines for dark sub- 
traction and flat fielding, spectral extraction, wavelength 
calibration with arc lamp lines, and the combining of 
multiple observations from a nodding sequence. The cor- 
rection f or atmospheric ab sorption lines was handled by 
xtellcor ([Vacca et al.ll2003l) . an IDL routine that divides 
the target spectrum by a calibrator. If the calibrator is 
an AO star, then its spectrum consists to a good approx- 
imation only of hydrogen lines and atmospheric absorp- 
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Table 1 

SpeX CTTS Sample 



Object Spectral Type 
(1) (2) 


log Afacc 

(3) 


Refs. 

(4) 


SpeX 
(5) 


HIRES 

(6) 


NIRSPEC 
(7) 


AA Tau . . 


K7 


-8.5 


5,2 


27 


30 


30 




K5 


-5.4 


3,3 


26,27 


30 


30,01 


BM And . . 


G8 


> -9 


6,1 


27 


30 


30 


BP Tau. . . 


K7 


-7.5 


5,2 


26 


30 


30 


CW Tau . . 


K3 


—6.0 


5,3 


26 


30,01 


30,01 


CY Tau . . 


Ml 


-8.1 


5,2 


27 


30 


30 


DF Tau^' . 


M2 


-6.9 


4,4 


26 


30 


30 


DG Tau A 


K7 


-5.7 


5,3 


26 


30 


30 


DK Tau A 


K7 


-7.4 


5,2 


26 


30 


30 


DL Tau. . . 


K7 


-6.7 


5,3 


27 


30 


30 


DO Tau . . 


MO 


-6.8 


5,2 


27 


01 


01 


DR Tau . . 


K7 


-5.1 


5,3 


27 


30 


30 


HN Tau A 


K5 


-8.6 


7,7 


27 


30 


30 


LkCa 8 . . . 


MO 


-9.1 


5,2 


27 


30 


30 


RW Aur A 


Kl 


-7.5 


7,7 


26 






UY Aur A 


MO 


-7.6 


4,4 


27 




01 



References. — (1) iQuenther fc Hessmai]| j993 : (2) IGHBC98I: (3) 
Hartigan ct al. 1995; (4) Hartigan & Kcnvon 2003; (5) IKenvorrTr^artmanril 
11995 : (6) .Rostopchinaia999. : (71 .White fc Ghez..20'5ll . 

Note. — Col. 2: Generally accurate to ±1 subclass; Col. 3: Logarithm of 
the mass accretion rate in Mq yr~^; Col. 4: References for the spectral type 
and mass accretion rate; Cols. 5—7: UT day of the month on which a spectrum 
was obtained in 2006: Nov 26, 27, or 30 or Dec 01. 

Unresolved binary (0.09" separation); properties are for the primary. 



tion lines. Division by such a calibrator, if observed at 
approximately the same airmass as the target, removes 
the target's atmospheric absorption lines. The hydrogen 
lines in the calibrator are removed by fitting a scaled 
model of Vega. After telluric correction, the merging of 
the six orders provided in the SXD mode and the removal 
of remaining bad pixels are accomplished with separate 
programs in the Spextool package. None of the targets 
showed extended line emission in their 2D spectra. 

The error in the flux calibration is dominated by dif- 
ferent slit losses for the target and the calibrator. We 
compared our approximate flux calibration to the J and 
K mag nitudes from 2MAS^]. Although we find that 

JspcX - •/2MASS = -0.11 ± 0.43 and i^SpoX - ^'^2MASS = 

—0.06 ± 0.41, the discrepancy is much less between our 
data and 2MASS in the J - K color, only -0.05 ± 0.20 
magnitudes, indicating that the magnitude differences 
are nearly color-independent, whether from variability 
or slit losses. Thus we are confident that despite uncer- 
tainties in the absolute flux calibration, the shape of the 
spectral energy distribution (SED) is preserved, and our 
continuum-normalized spectra fairly represent the wave- 
length dependence of the observed spectrum. 

2.2. Echelle Spectra: HIRES + NIRSPEC 

Beginning three nights after the SpeX run, high- 
resolution red optical and F-band echelle spectra of 14 
of the 16 CTTS from the SpeX sample were obtained 
simultaneously with HIRES on Keck I and NIRSPEC on 
Keck II, respectively. A 15th star from the sample was 
observed with NIRSPEC only. Stars with high-resolution 
spectra are noted in Table [TJ 

The Two Micron All Sky Survey (2MASS) is a joint project of 
the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the National 
Science Foundation. 



HIRES was used with the red collimator and the C5 
decker (1.15" x 7.0"), which has a projected slit width 
of 4 pixels and a spectral resolution of i? = 33, 000 
{AV = 8.8 km s~^). With the cross-disperser set to 
approximately 0.884° and the echelle angle at 0.0°, we 
have nearly complete spectral coverage from about 0.48 
to 0.92 fim in 38 orders (17 on the blue chip, 12 on the 
green chip, and 7 on the red chip). The red, green, and 
blue detectors were used in low-gain mode, resulting in 
readout noise levels of 2.8, 3.1, and 3.1 e~^, respectively. 
Internal quartz lamps were used for flat fielding, and 
ThAr lamp spectra were used for wavelength calibra- 
tion. The HIRES data were reduced with the MAuna 
Kea Echelle Extraction reduction script (MAKEE) writ- 
ten by Tom Barlow. Scott Dahm was the observer on 
Keck I during our HIRES run. 

NIRSPEC was used with the Nl fiher {Y band), which 
covers the range 0.95 to 1.12 /im in 14 spectral orders at 
a resolution R = 25, 000 {AV = 12 km s"^). Data reduc- 
tion, including wavelength calibration and spatial recti- 
fication, extraction of one-dimensional spectra from the 
images, and removal of telluric emission and absorption 
features, is performed with the IDL package REDSPEC 
by S. S. Kim , L. Pr ato, and I. McLean, as discussed in 
iFischer et afl (|200l . 

3. DERIVING THE EXCESS EMISSION SPECTRA 

We first determine the shape and intensity of the spec- 
trum of excess (non-photospheric) continuum emission 
between 0.48 and 2.4 /im using both echelle and SpeX 
spectra. We follow different procedures for the echelle 
and SpeX data and merge the results to give the SED of 
the excess emission over our full wavelength range. The 
starting point for both data sets is to measure individ- 
ual line veilings rx, defined as the ratio of excess flux 
E\ to photospheric flux P\ at a particular wavelength. 
Thereafter, we follow different routes to find the contin- 
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uous veiling spectrum V\ (a function describing r\ over 
a continuous span of wavelength), and the spectrum of 
the excess continuum E\. In both cases the excess con- 
tinuum is expressed in units of the photospheric flux at 
0.8 /im. Each step is described more fully in the following 
subsections. 

3.1. Echelle Spectra 

The determination of E\ from the echelle spectra fol- 
lows directly from the veiling in each order where a reli- 
able measurement can be made. To determine th e echelle 
line veilings r^, we follow iHartigan et al.l ()1991f) . match- 
ing each echelle order of an accreting star to that of a 
non-accreting standard of comparable spectral type to 
which has been added a constant excess continuum that 
reproduces the depth of the photospheric features in the 
accreting star. For this procedure, we use a modified 
version of an interactive IDL routine kindly provided by 
Russel White to accomplish the requisite velocity shift 
of the standard via cross-correlation and the requisite 
rotational broadening of the standard before determin- 
ing the required level of excess continuum. A polynomial 
is then fit to the individual r\ with wavelength, yielding 
V\ , and the excess emission E\ is the product of V\ and 
a temperature-matched photospheric template from the 
Pickles library ((Pickles 1998 ). 

Spectral templates for the line-veiling measurements 
include the w eak T Tauri star (WTTS) V 819 Tau (spec- 
tral type K7; iKenvon fc HartmannI 119951 ) and a grid of 
dwarf standards obtained with the same instrumental 
setup. The measurements of rx in the echelle spectra 
(HIRES and NIRSPEC) are shown in Figure [H along 
with the polynomial fits V\ to characterize the behavior 
with wavelength. In the optical spectra, veiling is mea- 
sured over a full order for those orders with strong pho- 
tospheric lines, no line emission, and no problems with 
corrections for telluric absorption, which are particularly 
acute at longer wavelengths. Although there are 38 opti- 
cal orders averaging about 100 A per order, the number 
of orders measured per star is considerably less than this 
and is not uniform among all the stars. Also, the order- 
to-order scatter is small in many cases, but in a few it 
is considerable. For the simultaneous NIRSPEC spectra, 
each of the 14 orders is about 120 A wide, and again, not 
all yield reliable veiling measurements. In the Y band, we 
find no appreciable changes in veiling with wavelength, 
so only a single value is plotted at 1.08 /^m in Figure [TJ 
Echelle veilings from both spectrographs are tabulated 
for nine wavelengths in Tabled where each HIRES value 
is an average of three adjacent orders (about 300 A) and 
the single NIRSPEC value is an average of several ad- 
jacent orders but is representative of the veiling over a 
range of about 1700 A. 

In the wavelength coverage of the echelle data, from 
0.48 to 1.1 /im, most stars show a monotonic decrease 
in veiling from the bluest wavelength down to 0.65 /im. 
Beyond that, some continue to decrease to longer wave- 
lengths, some fiatten out to a constant veiling approach- 
ing the 1 /im region, and two (DC Tau A and a second 
observation of CW Tau, shown in gray) rise toward 1 /im. 
The two stars with the lowest veilings, AA Tau and 
BM And, show no convincing wavelength dependence to 
the veiling through the optical region, and none of the 



stars show convincing differences in veiling through the 
Y band. The resulting spectra for the excess emission 
E\ will be shown in Section 3.3 along with the results 
from the SpeX data. 

3.2. SpeX Spectra 
3.2.1. SpeX Line Veilings 

For the lower-resolution SpeX spectra, line veilings r\ 
provide the basis for disentangling the effects of extinc- 
tion and veiling on the relatively flux-calibrated spectra, 
directly yielding the spectrum of the excess continuum 
emission Ex and the continuous veiling function Vx- 

The procedure to determine line veilings from the SpeX 
data is similar to that for the echelle spectra, although 
the process is hampered by the order-of-magnitude lower 
spectral resolution. For spectral templates we have 
dwarfs and gi ants in the IRTF SpeX spectral library 
(jRavner et al. 1 12009) plus our own SpeX spectra of three 
WTTS: HBC 407 (G8, observed on 2009 Dec 5), V819 
Tau (K7, observed on 2006 Nov 26), and LkCa 14 (MO, 
observed on 2009 Dec 30). We identified ten spectral seg- 
ments from 200 to 500 A wide that cover the strongest 
photospheric features, illustrated in Figure [2] for V819 
Tau. At SpeX resolution many of the features are un- 
resolved blends, but the strongest contributors to each 
region are identified in the figure. We found the WTTS 
provided better templates for line veilings than dwarf 
templates from the IRTF library as judged from a smaller 
scatter in rx with wavelength. We thus used HBC 407 
for BM And and V819 Tau for the rest of the sample. 

The veilings measured for each of the ten wavelength 
regions are listed in Table [3] For some of the high- veiling 
CTTS, either line or continuum emission made it impos- 
sible to determine veilings in all wavelength segments. 
The most severe cases were AS 353A and RW Aur A, 
where photospheric features could be discerned in only a 
few of the longer- wavelength regions near 2 /tm. In the 
spectrum of DC Tau A, the wavelength regions short- 
ward of 1 fim contain some photospheric features that 
are in emission. Although the specific lines used for veil- 
ing measurements display absorption, the adjacent line 
emission likely compromises the measurements in these 
regions. 

The presentation of the combined veilings from the 
echelle and SpeX data is in § 3.5, after Vx has been de- 
rived from the excess emission. Since this process de- 
pends on the low-resolution SpeX veilings, we need to 
assess their accuracy. We do this by comparing them to 
the echelle veilings in the region of spectral overlap. Dif- 
ferences could arise either from error in the SpeX veilings 
or from veiling variations in the 3-4 day interval between 
the two sets of observations. As a check on variability, 
we can also compare emission-line equivalent widths in 
the region of spectral overlap. 

These comparisons are made in Figure [31 The veiling 
comparisons are shown in the two upper panels, between 
SpeX and HIRES at 0.85 /im on the left and between 
SpeX and NIRSPEC at 1.05 /im on the right. The line 
equivalent widths are compared in the two lower pan- 
els, between SpeX and HIRES for Ca II A8500 on the 
left and between SpeX and NIRSPEC for Pa7 on the 
right. These equivalent widths, along with those for Br7 
from the SpeX data, are listed in Table H] (The num- 
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Figure 1. Echelle veilings for tiie 14 CTTS observed witli HIRES and NIRSPEC. Measured values are shown with - 
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ber of stars in each panel is less than the full 16 because 
two stars were not observed with HIRES, one was not 
observed with NIRSPEC, and shorter- wavelength SpeX 
veilings could not be measured in two stars.) 

In general we find good agreement between the veil- 
ings measured with the echelle and SpeX spectra taken 
a few days apart, although the correspondence is better 
between HIRES and SpeX than between NIRSPEC and 
SpeX for the majority of stars. Only three stars, all with 
high veilings, show significant veiling differences, up to a 
factor of 2, at these two wavelengths. One of them, DL 
Tau, shows higher veiling with SpeX than in the echelle 
spectra at both 0.85 and 1.05 /im, while the other two 
show discrepancies at only one wavelength. (DG Tau 
A has higher veiling at 0.85 /im with SpeX than with 
HIRES, and HN Tau A has higher veiling at 1.05 /xm 
with NIRSPEC than with SpeX.) As will be apparent 
when we introduce the full spectrum of the veiling Vx in 
§ 3.5, the excess emission of DL Tau clearly varied over 



a three-day interval, but for the other two stars the dis- 
crepant points are anomalously high compared to echelle 
and SpeX veilings at other wavelengths. We thus at- 
tribute them to error, possibly from line emission filling 
in photospheric features, and consider them to be upper 
limits. The line equivalent widths at the two epochs, 
which are not as sensitive to the difference in spectral 
resolution, are again similar for most stars, although 
Ca II A8500 is more variable on a three-day timescale 
than Pa7. We conclude that the lower-resolution SpeX 
line veilings in Tabic [3] are generally reliable, and we use 
them to derive the spectra of the excess continuum emis- 
sion in the next section. 

3.2.2. Method for Finding the Spectrum of the Excess 
Continuum from SpeX Spectra 

For the relatively flux-calibrated SpeX data, we use 
line veilings to extract the broad SED of the excess con- 
tinuum emission. We start with the assumption that 
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Figure 2. The ten spectral regions with the strongest photospheric features used to measure veilings in the SpeX data, shown for the 
WTTS V819 Tau. The dominant contributors to the absorption line spectrum in each region are indicated according to IRavner et al.l 

iHoog). 
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Table 2 

HIRES + NIRSPEC Average Line Veilings 



Object UT Date r4883 rsigs r5485 r5942 reses '"6725 ''7427 '"8102 rs695 riosoo 



AA Tau . . 


061130 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.0 


0.1 


0.1 


AS 353A . 


061130 


2.7 


3.3 


2.9 


2.4 


2.0 


2.0 


2.2 


2.3 


2.0 


1.7 


AS 353A . 


061201 




















1.4 


BM And . 


061130 


0.3 


0.3 


0.2 


0.2 


0.3 


0.2 


0.2 


0.2 




0.5 


BP Tau . . 


061130 


1.0 


1.0 




0.8 


0.6 


0.6 


0.6 


0.4 


0.4 


0.2 


CW Tau. . 


061130 


0.6 


0.9 


0.8 


0.7 


0.7 


0.8 


0.8 


0.7 


0.8 


1.7 


CW Tau. . 


061201 


1.0 


1.2 


1.1 


0.8 


0.9 


0.9 


0.9 


1.0 


1.2 


0.7 


CY Tau . . 


061130 


0.7 


1.0 


0.8 


0.5 


0.7 


0.6 


0.2 


0.0 


0.2 


0.2 


DP Tau . . 


061130 


2.0 


3.0 


1.7 


1.1 


1.7 


1.5 


0.9 


0.4 


0.6 


0.0 


DG Tau A 


061130 


1.7 


2.0 


1.4 


1.0 


1.0 


0.9 


0.8 


0.9 


1.1 


0.5 


DK Tau A 


061130 


1.0 


1.0 


0.7 


0.7 


0.7 


0.6 


0.5 


0.6 


0.4 


0.4 


DL Tau . . 


061130 


2.0 


2.0 


1.3 


1.0 


1.2 


1.3 


1.1 


1.4 


1.0 


1.0 


DO Tau . . 


061201 


2.5 


2.3 


2.0 


2.0 


1.2 


1.1 


0.8 


0.8 


0.6 


0.3 


DR Tau . . 


061130 


6.3 


7.3 


7.0 


6.0 


5.8 


5.0 


3.5 


3.0 


2.9 


3.5 


HN Tau A 


061130 


3.5 


3.3 


2.7 


1.4 


1.8 


1.6 


1.5 


1.0 


1.4 


(1.9) 


LkCa 8 . . . 


061130 


0.6 


0.6 


0.5 


0.4 


0.4 


0.2 


0.2 


0.1 


0.1 


0.1 


UY Aur A 


061201 




















0.3 



Note. — Parentheses indicate an upper limit (see § 3.2.1). 



Table 3 

SpeX Line Veilings 



Object 


UT Date 


ro.S2 


''0.91 


''0.97 


ri.os 


ri.is 


»'1.31 


'■1.98 


'■2.11 


^■2.20 


?'2.26 


AA Tau . . 


061127 


0.1 


0.25 


0.2 


0.2 


0.2 


0.2 


0.3 


0.3 


0.3 


0.5 


AS 353A . 


061126/27=" 
















3.8 






BM And. . 


061127 


0.3 


0.0 


0.2 


0.4 


0.5 


0.5 


1.4 


1.7 


1.8 


2.0 


BP Tau. . . 


061126 


0.3 


0.4 


0.4 


0.5 


0.5 


0.5 


0.8 


0.9 


1.0 


1.0 


CW Tau . . 


061126 


1.3 




1.3 


0.5 


0.9 


1.6 


4.5 


4.4 


8.1 


5.7 


CY Tau . . 


061127 


0.15 


0.2 


0.2 


0.2 


0.4 


0.3 


0.3 


0.5 


0.5 


0.6 


DP Tau . . 


061126 


0.5 


0.3 


0.4 


0.5 


0.8 


0.7 


0.7 


1.1 


1.1 


1.1 


DG Tau A 


061126 


(2.3) 




1.3 


0.9 


0.8 


1.5 


1.7 


1.2 


1.2 


2.1 


DK Tau A 


061126 


0.6 


0.6 


0.5 


0.5 


0.7 


0.6 


1.5 


1.9 


2.0 


2.0 


DLTau... 


061127 


1.8 




1.7 


1.8 


2.9 


1.7 


2.5 


2.3 


3.0 


3.1 


DO Tau . . 


061127 


0.7 




0.6 


0.7 


1.1 


1.3 


2.2 


2.6 


3.1 


2.9 


DR Tau . . 


061127 


3.0 




2.5 


3.5 




6.0 




9.0 


10 


10 


HN Tau A 


061127 


1.1 


1.0 


0.8 


1.0 


1.4 


1.5 


3.2 


2.4 


5.1 


3.2 


LkCa 8... 


061127 


0.0 




0.1 


0.1 


0.2 


0.2 


0.5 


0.6 


0.6 


0.7 


RW Aur A 


061126 
















3.9 




6.7 


UY Aur A 


()()1127 


O.r, 






O.T) 


0,8 


0..S 


l.i 


1,1 


1.9 


1.7 



Note. — Parentheses indicate an upper limit (see § 3.2.1). 
^ Two nearly identical spectra of 2006 Nov 26 and 27 were averaged to increase the signal-to-noise 
ratio. 
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Figure 3. Comparison of veilings and emission-line equivalent widths from SpeX and HIRES+NIRSPEC, separated by 3-4 days, in the 
region of overlap. Top left: Veiling near 0.85 /im in the HIRES versus the SpeX data. Top right: Veiling near 1.05 fivci in the NIRSPEC 
versus the SpeX data. Bottom left: Equivalent width of the Ca II A8500 line in the HIRES versus the SpeX data. Bottom right: Equivalent 
width of the Pa7 line in the NIRSPEC versus the SpeX data. 



the spectrum 0\ of a CTTS observed with SpeX can be 
described as the reddened sum of the flux from the pho- 
tosphere Pa and the flux from a continuum excess Ex, 
such that 



Ox = {Px+Ex)lO 



PA(l + rA)10 



(1) 



To find Ex we apply the method of IGHBC981 used 
to derive excess spectra from medium-resolution optical 
spectrophotometry of CTTS. The intrinsic CTTS pho- 
tosphere Pa will differ from that of a spectral template 
Pa of identical temperature with extinction Ax^t by a 
wavelength-independent scaling factor C that reflects the 
imprecision in the absolute flux calibration (§[23J as well 
as the different distances and radii of the photosphcric 
standards and the program stars: 



Pa = CTaIO"-*-^^-^. 



(2) 



The extinction to the program star is then determined 
by evaluating equation ([ij at the set of wavelengths 
where individual line veilings rx have been measured. 



Via substitution, equation ([T]) can be rewritten as 

Ox = CTa (1 + rx) 10-"-4(^a,o-a^,t) 

= CTa (1 + rx) lO-"-4feA(^v,o-Av,T)^ (3) 

where the extinction law kx — Ax/ Ay is assumed to be 
th e same for b o th obj ects. Here we use the extinction law 
of lFitzpatricl3 (|l999l i with Ry = 3.1, as represented in 
the routine fm_unred.pro in the IDL Astronomy Library^ 
Reformatting by rearranging terms, taking the loga- 
rithm of both sides, and multiplying by 2.5 yields 

2.5 log [Pa (1 + rx) /Ox] = fcA {Av.o ~ Av,t) - 2.5 log C, 

which is identical in form to equation (5) of IGHBC981 
This equation is linear in fcA, such that a line described 
by it has slope Ay^o — Av,t and intercept — 2.51ogC. 
Thus both of these are readily determined from a linear 
fit to2.51og[PA(l + rA)/OA] 

vevsus Ax/Ay. Once Ay o, 
the extinction to the CTTS, is found, is it straightforward 
to recover the spectrum of the excess emission from the 

* http://idlastro.gsfc.nasa.gov/ 
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Table 4 

Emission-Line Equivalent Widths 





SpeX 


HIRES 


SpeX 


NIRSPEC 


SpeX 




Ca II A8500 


Ca II A8500 


Pa7 


Pa7 


Br-v 














AA Tau . . 


0.0 


-0.4 


0.9 


0.3 


1.6 


AS SSSA'' 


49.6 


52.9 


14.8 


15.5 


17.9 


AS 353A . 


53.2 




18.4 


14.7 


21.3 


BM And . . 


-0.9 


-0.4 


-0.8 


-0.1 


0.9 


BP Tau. . . 


3.5 


3.0 


6.3 


4.5 


4.7 


CW Tau'' . 


9.9 


10.8 


5.2 


4.3 


3.3 


CW Tau . . 




11.0 




5.7 




CY Tau . . 


0.0 


3.3 


1.2 


4.3 


1.1 


DF Tau . . 


2.1 


3.1 


4.2 


3.3 


3.3 


DG Tau A 


39.3 


39.9 


9.3 


8.5 


7.4 


DK Tau A 


2.8 


1.6 


4.2 


2.5 


2.2 


DL Tau. . . 


33.3 


49.9 


16.4 


16.4 


12.0 


DO Tau . . 


12.6 


27.4 


6.7 


8.8 


2.5 


DR Tau . . 


39.4 


25.5 


18.4 


11.2 


8.6 


HN Tau A 


26.8 


46.9 


7.7 


8.3 


4.3 


LkCa 8 . . . 


0.0 


0.2 


1.2 


0.6 


1.5 


RW Aur A 


69.6 




13.7 




10.2 


UY Aur A 


2.0 




3.9 


1.4 


2.0 



^ AS 353 A and CW Tau were observed twice by at least one instrument. 
The data from the epochs marked with this symbol are used in all plots. 



reddened and veiled Ox- 

The key steps in the above process of extracting the 
SED of the excess emission are illustrated in Figure |4] for 
BP Tau. The upper panel shows the observed spectra of 
both BP Tau, X0\, and the K7 V photospheric template 
HD 237903 from the IRTF library, ATa, plus the ten mea- 
sured line veilings rx . The spectra are plotted in units of 
XFx scaled to unity at 0.8 fj,m, although the spectrum of 
BP Tau is shifted upward by 0.25 for clarity. The second 
panel plots 2.5 log [Tx {I + rx) /Ox] against kx at each of 
the ten wavelengths for which the veiling has been mea- 
sured. As indicated in equation (|4]), this should be a lin- 
ear relation, with a slope equal to A vq — Avt- Since th e 
template HD 237903 has Av,t = (|Ravner et al.ll2009D . 
the slope of the best-fit line is the extinction toward the 
object Av.o- (Note that, because kx increases from right 
to left, a line that increases from right to left represents 
a positive Ay^o-) Once the extinction is found, the ob- 
served spectrum XOx can be dereddened, as shown in the 
third panel. If the template spectrum is then multiplied 
by C (the normalization constant corresponding to the 
intercept of the best-fit line) , it can be subtracted from 
the dereddened object, yielding the excess spectrum XEx, 
also shown in panel 3. The template spectrum in panel 
3 is scaled such that XFx = 1 at 0.8 /im, and the other 
two spectra are shown with the correct scaling relative to 
the template. We will express the excess XEx in units of 
the photospheric continuum at 0.8 /im in all subsequent 
figures and tables. 

Finally, the continuous veiling spectrum Vx, shown in 
panel 4, is found by dividing the spectrum of the excess 
emission Ex by the scaled photospheric template Px ■ We 
compare it to the measured line veilings rx as a consis- 
tency check. If the line veilings lie close to the derived 
veiling spectrum, as seen for BP Tau in Figure IH then 
the template is a good match to the CTTS photosphere, 
and Ay is well determined. 

The normalization constant C is related to the ra- 
dius of the T Tauri star if the absolute flux calibration 
is accurate. Provided the spectra of the template and 



the dereddened, deveiled CTTS have the same wave- 
length dependence, their fiux ratio can be written C = 
e(i?cTTs/'^CTTs)^/(-Rtcmp/rftemp)^, where e accounts for 
errors in the absolute flux calibration. Solving for the 

CTTS radius, i?CTTS = (C/e)l/2 (dcTXs/rftomp) i?tcmp. 

With BP Tau as an example and HD 237903 as 
its template, C = 6.4 2 x 10~^; dtcmp = 12.9 pc 
(jCould fc Chanainll2004D : rfc TTS = 140 p c, a commonly 
accepted distance to Taurus (iKenvon et al. 2008); and 
R^^^p = 0.67 Rq (| Johnson fc Wrightlll983i) . giving a ra- 
dius of 1.85 i?0±2O% for BP Tau if the error is dominated 
by a 40% uncertainty in the absolute flux calibration 
(§ 2.1). This compares favorably with other estimate s 
for the radius of this star: 1.9 Rq (iHartiean et al."1995f). 
L99 i?0 (GHBC98i) . and 1.95 Rr:. (IJohns-KruU et all 
|1999( ). indicating that either the SpeX absolute flux cali- 
brations for BP Tau and its telluric standard are reason- 
ably accurate, or that any inaccuracies cancel when the 
CTTS is divided by its telluric standard. 

Before showing the results of this procedure, summa- 
rized in Figure m to find Ay, the spectrum of the excess 
emission Ex , and the spectrum of the veiling Vx for each 
star in the SpeX data set, we demonstrate how the out- 
come depends on the properties of the chosen extinction 
law and spectral templates. 

3.2.3. Sources of Uncertainty in Determination of Av and 
Ex from SpeX Spectra 

The extraction of the excess emission spectrum from 
our SpeX data depends on the assumed extinction law 
and the spectrum of the templ ate. We experime nted 
with several e xtinct ion laws from iFitzpatrickl ()1999[ ) and 
iCalvet et al.l ([20041) thought to be more appropriate for 
dark clouds than the standard ISM extinction with Ry — 
3.1. We compare in Figure [5] the effect on the derived 
Ay and ij^A between extinctio n laws with Ry = 5.0 and 
Ry = 3.1 (|Fitzpatricklll999f) . The larger Ry yields a 
small decrease in Ay of 0.2 mag (1.55 instead of 1.75, 
a 10% effect) and an excess spectrum that differs from 
the Ry = 3.1 result by an average of 3% (up to 6% at 
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Figure 4. Method for deriving Ay, the excess eiiiissioii spectrum Ex, and the continuum veiling Vx with SpcX data. Top panel: Observed 
spectra of the object BP Tau (black and offset for clarity), the template K7 dwarf HD 237903 (gray), and the SpeX line veilings rx for BP 
Tau (+). Second panel: Veiling-corrected logarithm of the template- to-object flux ratio at each wavelength with a measured I'x, plotted 
against the ratio of Ax to Ay for the adopted reddening law (see eqn. 4). The slope of the best-fit line is Ay o — Ayj-. Third panel: 
Normalized, reddening-corrected spectra of BP Tau (black), the K7 V template (gray), and their difference (black), i.e., the excess spectrum 
XEx. Bottom panel: The continuous veiling (solid curve), resulting from dividing Ex by the template, and the line veilings rx (+)■ 
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any one wavelength). The effects are small in the SpeX 
SXD regime, since the extinction laws have nearly the 
same ratios of A\/Av from 1.2 to 2.5 /im (0.27 at J, 
0.17 at H, and 0.11 at K) and differ only modestly from 
0.8 to 1.2 fim. We thus adopt the standard interstellar 
extinction law (Ry = 3.1) with the consequence that Ay 
may be overestimated by a few tenths of a magnitude; the 
derived spectrum of the excess emission is not affected 
significantly. 

The derivation of Ay and E\ are also affected by the 
choice of a spectral template. Ideally the photospheric 
template Tx would be matched in both temperature and 
surface gravity to the CTTS Ox in equation 4, and its 
extinction would be perfectly known. Thus a good grid 
of WTTS with well defined extinctions would be prefer- 
able templates, but we acquired only three SpeX spectra 
of WTTS, with spectral types G8, K7, and MO. An al- 
ternative is to use the IRTF library, where we can select 
unreddened dwarf templates to match the spectral type 
of each program star. 

In Figure|6]we illustrate the consequences in the deriva- 
tion oi Ay, E\, and Vx for BP Tau that result from differ- 
ent choices for the spectral template. In the top panel, we 
show the effect of errors in the temperature classification 
of BP Tau (nominally K7) by plotting its XEx with three 
dwarf spectra from the IRTF library as templates: a K7 
and the next hottest (K5) and next coolest (MO) dwarfs 
available in the library. While spectral types of CTTS 
are known to ±1 subclass (see references in Table [T]) and 
those of the templa tes are known to ±0.65 subclasses 
(|Ravner et al.ll2009D . the fi gure shows that a mismatch 
as large as two subclasses (K5 versus K7) yields an error 
in the excess spectrum of only about 10% at the most 
discrepant wavelengths. Although the Ay derived vary 
substantially (2.17 when the K5 dwarf is used, 1.75 when 
the K7 is used, and 1.48 when the MO is used) this uncer- 
tainty does not impact our results on the emission excess 
because we use a reddening- independent approach to find 
E\ from our optical spectra. 

In the middle panel we show XEx from three templates 
with the same temperature and different surface gravi- 
ties: (1) the K7 V from the IRTF library, (2) a Kurucz 
model with log.9 = 5.0 and T = 4000 K, and (3) a Ku- 
rucz model with log 17 — 4.0 and T — 4000 K. Since the 
Kurucz model^ have lower spectral resolution than the 
SpeX data, the continuum shape from each model was 
combined with the lines from the IRTF standard so pho- 
tospheric features would subtract cleanly. The compari- 
son shows that similar excesses are derived from both the 
IRTF template and the logg = 5.0 model. The excess 
derived from the log 5 = 4.0 model is slightly smaller in 
the J and K bands and more markedly so in the / and H 
bands. In particular, the derived excess from the lower 
surface-gravity model is flat between 0.8 and 0.85 /im and 
through the H band, while for the higher surface-gravity 
templates the spectrum rises shortward of 0.85 fim and 
shows a small hump at H. This comparison suggests 
that if the Kurucz log g = 4.0 models are more appropri- 
ate templates for CTTS, then our use of dwarf standards 
will introduce a small error into the derivation of Ex. 

In the bottom panel of Figure [BJ we compare XEx for 
BP Tau using two WTTS templates with spectral types 

^ ftp:/ /ftp. stsci.edu/cdbs/grid/ck04models/ 



K7 and MO to the result from the IRTF K7 V stan- 
dard. One WTTS, LkCa 14 (MO), has zero extinction 
and yields a XEx for BP Tau that shows a much better 
match to the excess found with the K7 V template than 
with the Kurucz log 5 = 4.0 model. The other WTTS, 
V819 Tau, of spectral type K7, does not have negligible 
extinction. If we adopt Ay — 2.6 for V819 Tau, which 
follows from dereddening it against either LkCa 14 or a 
K7 V (see next paragraph), the derived excess for BP 
Tau is again comparable to that from LkCa 14 and the 
IRTF K7 V dwarf, although it is flatter in the /, H, and 
K bands, and it is not a good match to the excess from 
the Kurucz \ogg = 4.0 model. 

The comparison shown in Figure [6] suggests that the 
IRTF library, offering a flne grid of dwarf spectral types 
with no extinction, provides a reasonable choice for our 
spectral templates. A potential consequence of this ap- 
proach is that the surface gravity of a dwarf, typically 
log 5 = 4.5, is somewhat higher than the log (7 = 4.0 
found for f our WTTS from hig h-resolution photospheric 
line ratios (jTakagi et al.l 120101 ) . but not out of line con- 
sidering the results in Figure [5] and the sizable scatter in- 
ferred for log g distributions of pre-main-seq u ence stars 
from colors and line ratios (Hille nbrandll2009( ). Surpris- 
ingly, however, no matter which spectral template we 
use, we derive an extinction to BP Tau that is si gnifl- 
cantly larger than the Ay of 0.5 found bv lGHBCOl . Our 
values for Ay range from 1.75 to 1.05 (see above and 
caption to Fig. |6|), all larger than the overestimate of a 
few tenths of a magnitude that would be expected if we 
adopted a dark-cloud extinction law. 

The large Ay we consistently find from comparing 
CTTS to spectral templates in the SpeX SXD spectral 
regime may be related to a peculiar beh avior i n WTTS 
first noted by GHBCOj and ICullbrin g eTall (|l998b). 
They found that one of the WTTS templates we use 
here, V819 Tau, has anomalous colors, becoming increas- 
ingly redder with wavelength than a main-sequence star 
of comparable spectral type. The observed color anoma- 
lies in V819 Tau, also found in a number of other WTTS 
but not in our other SpeX template LkCa 14, lead to 
the curious result that Ay derived from colors increases 
with wavelength, with Ay = 0.5 from B~V to Ay — 1.8 
from V - J. Indeed, the Ay = 2.6 we found for V819 
Tau when dereddening its SpeX spectrum against either 
the WTTS LkCa 14 or a K7 V template (both with no 
reddening) continues this trend. GuUbring et al. consid- 
ered possible causes of anomalous colors and favored the 
(untested) explanation that very large cool spots might 
produce the observed effect. That V819 Tau has an un- 
recognized emission excess in the spectral region of our 
study is worth considering, since it was recently found 
with Spitzer to h ave a weak IR ex cess at wavelengths 
exceeding 10 /xm (|Furlan et al.ll2009( ). 

As will be shown the next section, the high Ay de- 
rived from SpeX SXD data is not unique to BP Tau; 
it is also the case for most of our CTTS. Fortunately, 
this does not have a significant effect on our derived ex- 
cesses, since in the SpeX wavelength region the extinc- 
tion correction is not large, and in the optical we use a 
reddening-independent method to recover the SED of the 
excess emission. We return to this topic in the discussion 
section. 

We conclude from the good agreement between the 
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Figure 5. Comparison of extinction laws witii different Ry in the derivation of the excess spectrum Ex of BP Tau. Top panel: Veiling- 
corrected logarithm of the template-to-object flux ratio at each wavelength with a measured rx, plotted against the ratio of Ax to Ay for 
each extinction law. The slope of the linear relation is Ayp — Ay j- (eqn. 4). Bottom panel: The corresponding excess spectrum XEx for 
each extinction law. Although the derived Ay is smaller for Ry = 5.0, the difference in the derived XEx is not significant at SpeX SXD 
wavelengths. 



XEx derived from the unreddened WTTS LkCa 14 and 
the unreddened K7 V shown in Figure|6]that using unred- 
dened dwarfs from the IRTF hbrary matched in spec- 
tral type to each CTTS is a reasonable choice for spec- 
tral templates, a practice also adopted bv lEspaillat et al.l 
dfOlO) in the near infrared. However, the shape of the 
derived XEx in the / and H bands may be slightly af- 
fected if there is a gravity mismatch between the CTTS 
and the template. 

3.2.4. Derivation of Ay from SpeX Spectra 

In this section we derive the extinctions to the CTTS 
program stars from their SpeX spectra and observed 
line veilings following the method demonstrated in Fig- 
ure m For the reasons outlined in the previous sec- 
tion, we use IRTF dwarf standards with no extinction 
as spectral templates, and we use a standard extinction 
law. The validity of these assumptions can be exam- 
ined in Figure [71 where we show the relation between 
2.5 log \Tx (1 + rx) /Ox] and kx (see eqn. 4 and panel 2 
in Fig. 13]), expected to be linear for each star. For most 
stars there is little scatter around the best-fit line, indi- 
cating that the SpeX line veilings are robust. The slope 
of each line corresponds to the difference in extinction be- 
tween the standard and the template, but since we have 
chosen IRTF standards with zero reddening, it gives Ay 
for each star, as identified in each panel and listed in Ta- 
ble [S] The procedure had to be modified for the two stars 
for which only one or two veiling measurements could be 
determined due to line emission filling in the veiled pho- 



tospheric features (AS 353A and RW Aur A). For these 
stars we used the published Ay to set the slope of the 
line, anchoring it through the one or two valid veiling 
measurements . 

One of our CTTS, AA Tau, is subject to periodic oc- 
cultation by its circumstellar disk, and thus it may have 
anomalous extinctio n if seen throu g h the disk. Following 
the prescription of iBouvier et al.l ()2007D , in which the 
period of AA Tau is reported to be 8.22 days and JD 
2453308 corresponds to phase 0.51, our Keck data were 
acquired at phase 0.21, and our IRTF data were acquired 
at phase 0.85. As photometric and spectroscopic diag- 
nostics of disk occultation appear only between phases 
0.3 and 0.8, it seems we are not sampling extinction 
through the disk. The low veiling we observe (0.1 in the 
optical and 0.3 at K) is consistent with the level found 
by Bouvier et al. outside of the occultation phase, further 
supporting this conclusion. 

We also show in Table [5] both the normalization con- 
stant C that scales the template to the dereddened, de- 
veiled CTTS and the corresponding stellar radius for 
each CTTS. As described in §3.3, C depends on the 
radius of and distance to both the CTTS and the 
template. We adopted a distance to the Taurus ob- 
j ects of 140 PC a distance to AS 353A of 200 pc 
([Rice et al.l 12006^ a distance to BM And of 440 pc 
(jAveni fc Hunted 1 19691) . d istances to temp lates other 
than H P 237903 IC ould fcChanamel [200l from Hip- 
parcos (jPerrvman et al.l 1199'^ . and template radii from 
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Figure 6. The excess spectrum \Ex of BP Tau (K7) for different choices of template. Top: XEx based on IRTF K5 V, K7 V, and MO V 
standards. The resulting Ay for BP Tau are 2.17, 1.75, and 1.48, respectively. Middle: XEx repeated for the same K7 V standard and 
for two Kurucz models (adjusted to represent SpeX resolution; see text) with T = 4000 K and logg = 5.0 and 4.0. The resulting Ay for 
BP Tau with these Kurucz models as templates are 1.58 and 1.05, respectively. Bottom: XEx repeated for the same K7 V standard and 
for two WTTS, LkCa 14 (MO) and V819 Tau (K7), with resulting Ay for BP Tau of 1.71 and 1.75, respectively The XEx are normalized 
such that the photospheric XFx = 1 at 0.8 fj,m. 



iJohnson fc Wrightl ([TOSl . The CTTS radii derived from 
C lie between 0.8 and 3.4 Rq with a mean of 1.8±0.7 Rq. 
As this is a typical CTTS radius, we infer that the ab- 
solute flux calibration of the SpeX data is reasonably 
accurate. 

The extinctions we derive from the SpeX data are typ- 
ically 0.5 to 1 magnitude higher than those from pub- 
lished optical data, as shown in Table [SJ As described in 
the previous section, a modest overestimate of Ay could 
come from using an extinction law that is not appropri- 
ate for dark clouds or from a spectral type mismatch, 
but these cannot account for the large differences found 
here. Although our results are not compromised for the 
reasons given above, the tendency to derive larger Ay 
in the SpeX regime for some WTTS and most CTTS is 
a subject for continued investigation, best addressed by 
a systematic campaign to acquire a dense grid of flux- 
calibrated WTTS spectra with wide wavelength cover- 
age. 



3.3. The Veiling and Emission Spectra from 
0.48 fim to 8.4 fJ-m 

In this section we present the spectra of the contin- 
uous veiling Vx and the continuous excess emission Ex 
obtained by merging the SpeX results with the HIRES 
and NIRSPEC results, thus covering the full wavelength 
range of our data, from 0.48 to 2.4 /im. 

The results for Vx are assembled in Figure 13 showing 
(1) the individual line veilings r\ from both the SpeX and 
echelle spectra, (2) the SpeX veiling curves Vx derived 
from the emission excesses, which include line emission, 
as in the fourth panel of Figure HI and (3) the echelle 
veiling curves Vx that are polynomial fits to the echelle 
veilings. In most cases there is reasonable agreement in 
the region of spectral overlap between the low- and high- 
resolution data, taken several days apart. The exception 
is DL Tau, where the SpeX V\ is about twice as high 
as the veiling from the echelle data; we attribute this 
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Table 5 

Results of Fits to SpeX Line Veilings 



Figure 7. Derivation of Ay (same as the seeond panel of Fig. \^ for each star in the sample. The slope of the best-fit line to the 
veiling-eorrected flux ratios ('+') is Ay, since in all cases the template has no extinction. The two stars with an asterisk by their Ay 
(AS 353A, RW Aur A) use published Ay for the slopes, since there are not enough rx points to derive them self-consistently. The spectral 
type of the IRTF dwarf template used for each star is given in parentheses. Stars are sorted by their 1 /im excess. 

are illustrated in Figure [9l plotted as XE\ in units of the 
photospheric flux at 0.8 /.tm. Over the range observed 
with SpeX, between 0.80 and 2.43 //m, the E\ c urves 
are generated with the procedure discussed in § 13.2.21 
and illustrated in the third panel of Figure SI which 
yields emission lines as well as the continuum shape. 
For the echelle data, between 0.48 and 1.08 fim, the 
XEx curves are the product of polynomial fits to the 
line veilings (Fig. [S]) and polynomial fits to the continua 
of temperature-matched dwarf standards from the Pick- 
les library (Pickles 1998): this method is independent of 
extinction. For comparison, we also show (plus signs) 
the point-by-point products of individual veiling mea- 
surements and the corresponding Pickles template at the 
same wavelength. This shows the scatter in the indi- 
vidual line-veiling measurements, including fluctuations 
that would result from measuring within broad molec- 
ular bands. In most cases the individual measurements 
show little scatter around the curve. However, for the 
few stars that show substantial scatter (e.g., CY Tau 
and DF Tau), the curves for XE\ may be less accurate. 
(The scatter for DF Tau may be due to its status as an 
unresolved binary.) In most stars there is good agree- 
ment between the excesses from the echelle and SpeX 
regimes, indicating that any systematic errors between 
the two approaches for finding E\ are not large. 

To facilitate evaluation of the excess emission among 
the sample objects, we compute a median excess E^x 
in each of seven wavelength intervals AA. The intervals. 
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difference to variability (see § 3.2). 
The resulting SEDs of the continuum excess emission 
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Figure 8. Line-veiling measurements rx compared to derived continuum veilings V\ over the wavelength range 0.48-2.4 fim with stars 
sorted by their 1 fim excess. Points for rx are (+) for HIRES and NIRSPEC and (x) for SpeX. Curves for Vx are black for HIRES and 
NIRSPEC (fits to individual rx) and dark gray for SpeX data (as in the bottom panel of Fig.|4]l. A second echelle spectrum for CW Tau 
is overplotted in light gray. 



defined in Table [HI roughly follow the passbands of stan- 
dard photometric filters V, R, I, Y, J, H, and K. The 
corresponding E/\\ in each wavelength interval are listed 
in Table [6] for each star, where Ey through Ej are from 
the echelle spectra, and Ey through Ek are from the 
SpeX spectra. As for XE\ in Figure El the median excess 
Eax in each wavelength interval is in units of the pho- 
tospheric flux at 0.8 /im. These median values will be 
useful in characterizing the lYJ excesses of CTTS that 
span several orders of magnitude in disk accretion rate. 

4. BEHAVIOR OF THE EXCESS EMISSION SPECTRA 

The CTTS in our sample show excess emission over the 
full range of our wavelength coverage. Although the me- 
dian excesses in the seven wavelength bands approximat- 
ing V, R, I, Y, J, H, and K oversimplify the full SEDs, 
a plot of the E/^\ sequence for all 16 stars in the sam- 
ple, shown in Figure [TOj is instructive. The overall shape 



of the excess emission is similar among the full sample, 
where, in units of the photospheric flux at 0.8 /im, the 
excess declines from V through Y io K while covering 
a range in magnitude at Y from three times the photo- 
spheric flux (DR Tau) to only a few tenths of the photo- 
spheric flux (AA Tau, LkCa 8). This can be shown quan- 
titatively by plotting the median excesses at 1 //m, Ey , 
against the median excesses at the shortest and longest 
wavelength intervals in our dataset, Ey and Ek- Com- 
paring the median excesses in these three bands is par- 
ticularly instructive, since Ey will include a contribution 
from the blue excess used to model disk accretion rates, 
Ek will include a contribution from the near-infrared ex- 
cess attributed to sublimating dust in the disk, and Ey 
will have minimal contributions from these two compo- 
nents. This comparison is shown in Figure [TTJ both for 
the excesses Ey, Ey, and Ek and for the correspond- 
ing continuum veilings VV, Vy, and Vk- The ratios of 
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Figure 8. cont. 

these excesses are remarkably consistent from star to star 
when plotted logarithmically. The median values of the 
excess ratios are Ey/Ey — 0.50 and Ey /Ek — 1-45, and 
the median values of the veiling ratios are Vy /Vv — 0.52 
and Vy /Vk = 0.35. The figure and columns 9 and 10 of 
Table ini show that the dispersions of the excess ratios for 
individual stars around these median ratios are not large; 
the interquartile ranges are 0.42 to 0.66 for Ey/Ey and 
1.25 to 2.20 for Ey/EK- 

Despite the correlations in Ey, Ey, and Ek shown in 
Figure [m the XE\ curves of individual stars in Figure IHl 
have different slopes in particular wavelength intervals 
that are masked when we take an average across each 
band. We illustrate this in Figure [12] for the J and K 
bands, since the spectra there are usually well described 
by a single slope. The slopes, included in Table HI are 
defined as ^aa = (A/E'a/ — KE\-) /\iE\., where \i and 
\f are, respectively, the shortest and longest wavelengths 
in the regions defined in Table [S] The slopes in both 



bands run roughly from —0.1 to 0.25 but show no cor- 
relation with each other or with the magnitude of Ey. 
The most notable outliers are three stars with large Ey 
and negative Sj: DL Tau, RW Aur A, and AS 353 A. 
Looking back at Figure [HI these three stars are distinc- 
tive, with XE\ rising steeply to shorter wavelengths from 
the J through the / band. In contrast, the other stars 
with large Ey are much flatter through the / band and 
have excesses that rise to longer wavelengths through the 
J band, with slopes from 0.1 (DR Tau, HN Tau A, and 
DG Tau A) to 0.2 (DO Tau and CW Tau). 

The two wavelength regions sensitive to surface grav- 
ity also show a range of behaviors in XEx (see § 3.4 and 
Fig. [6]). In the H band, a hump in XE\ becomes more 
prominent in stars with small excesses, as might be ex- 
pected from a surface gravity mismatch. However, since 
a rise to shorter wavelengths between 0.80 and 0.85 ^m, 
which would also be an expected from a surface gravity 
mismatch, is not seen in these stars, the structure in the 
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Figure 9. Emission excess curves XEx over the wavelength range 0.48-2.4 ura with stars sorted by their 1 /im excess. The SpeX excesses, 
which include emission lines, are calculated as in the third panel of Fig. g] For the HIRES and NIRSPEC excesses, the solid curve is 
the product of the fit to the line veilings rx and the fit to the continuum of the photospheric template, and plus signs are point-by-point 
products of each individual line veiling and the flux in the template at that same wavelength. Excesses are in units of the photospheric 
\Fx at 0.8 /im. A second echelle spectrum for CW Tau is overplotted in gray. 



H band may be inherent in the excess. 

Only a few stars have excess spectra that rise to shorter 
wavelengths between 0.80 and 0.85 /im (e.g., DL Tau). 
While this behavior may be related to the Paschen jmxip 
at 0.8207 /xm, no star shows a distinct Paschen dis- 
continuity in Ex . Figure [T3] plots an index measuring 
the Paschen jump (the ratio of the excess emission at 
0.815 /im to that at 0.915 /im) against Ey- In contrast 
to the Balmer jump, whic h is more pronounced for star s 
with weaker blue excess (|Herczeg fc Hillenbrandl [20081 ) . 
there is no correlation between the Paschen jump and 
the lYJ excess. 

We next compare the average excess continuum emis- 
sion in the three wavelength intervals V, Y, and K to 
the fluxes in strong emission lines that have been shown 
to correlate w ith non-simultaneously determined disk ac- 
cretion rates (jMuzeroUe et al.|[l998al lH). In Figure [HI we 



compare the excess continuum fluxes Ey, Ey, and Ek 
to veiling-corrected equivalent widths for the three emis- 
sion lines Ca II A8500, Pa7, and Br7. The directly mea- 
sured equivalent widths (see Table [4]) are multiplied by 
the factor 1 + r\ to make them relative to the photo- 
spheric flux, thus facilitating comparisons of stars with 
different levels of veiling. Although fluxes in each of the 
three lines correlate with the median excesses, the corre- 
lations show much more scatter for Ek than for Ey or 
Ey. The linear correlation coefficients, listed in the up- 
per left corner of each panel, are lowest for Ca II A8500 
in all three wavelength intervals, but this arises from a 
break in the dependence of the Ca II A8500 flux on the 
excess strength, which is weak at low excess and stronger 
at high excess, such that a single linear fit does not rep- 
resent the situation adequately. As will be discussed in 
a forthcoming paper examining the full complement of 
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Figure 9. cont. 



emission lines in our combined HIRES, NIRSPEC, and 
SpeX data set, this break in the Ca II versus Eax re- 
lation corresponds to a transition from the line being 
dominated by a narrow component to being dominated 
by a broad component. In contrast, the correspondence 
of the h ydrogen lines, which always show a broad com- 
ponent fEdwa rds et al.|[2006f) , to Ey or Ey is essentially 
linear, although the better correlation with Ey suggests 
that the Y band excess may be more closely affiliated 
with the region responsible for the hydrogen line emis- 
sion. 

In summary, we find that the levels of the excess emis- 
sion in three wavelength intervals that correspond to the 
beginning, middle, and end of our spectral coverage {Ey, 
Ey, and Ek) are well correlated with each other and 
with emission line fluxes recognized as accretion diag- 
nostics. Furthermore, while the shape of the excess over 
the full wavelength range is not strictly uniform from 
star to star (e.g., slopes vary over restricted wavelength 



intervals), it is roughly similar from star to star, making 
it unlikely that random effects such as large cool spots 
or undetected companions are major influences on the 
derivation of the excess. 

5. MODELING THE lYJ EXCESS EMISSION 

In this section we present an initial characterization of 
the temperature, size, and luminosity of the region(s) 
responsible for the excess emission between 0.48 and 
2.4 /im. According to previous investigations of the op- 
tical/UV and near-infrared excesses, the excess in this 
wavelength regime will have a contribution at its short- 
est wa velengths from hot accretion spots on the star 
(|CG98ll and contributions at its longest wavelengths 
from the dust sublimation zone in the disk |MCHD03). 
Based on these studies, we adopt Thot = 8000 K and 
?cooi = 1400 K for these two regions. The magnitudes of 
their contributions will then depend on the area of each 
region projected along the line of sight. For both regions, 
we express this area as a filling factor / relative to the 



19 



Table 6 

Median Excess Fluxes, Ratios, and Slopes 
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Note. — E^x ^■re the median excesses -Fa, total ~ ^A, photosphere over the wavelength ranges whose extents are 
indicated in microns, scaled so that Fx photosphere = 1 at 0.8 /im. Sa\ are slopes within the indicated bands; see 
§4. 

^ For DL Tau, Ey, Er, and Ej are multiplied by 1.7 to correct for a shift between optical and IR data attributed 
to variability. 
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Figure 10. Comparison of the median excesses i?AA for all stars in our sample in wavelength intervals that approximate the bandpasses 
V, R, I, Y, J, H, and K as defined in Table|5] This simplified plot of the emission excesses from Fig. |§]allows a comparison across the full 
sample of CTTS. All excesses are in units of the photospheric flux at 0.8 /tm. 
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Figure 11. Relations among excess emission at V, Y, and K (Table|6ll expressed both as veiling Vaa and emission -Eaa- Top: Vy versus 
Vv [left) and Vk {right). Bottom: -Eaa for the same pairs of wavelengths in units of the star's photospheric flux at 0.8 lira. Dotted lines 
correspond to the median ratios labeled in each panel. To account for variability, the V excess of DL Tau has been increased by the ratio 
of its SpeX to its echelle excess in the region where they overlap. 
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Figure 12. Slopes of \Ex in the J band {left) and the K band (right) versus the normalized flux of the continuum excess in Y. The slope 
across a band AA is defined as Saa = i^f^Xf ~ ^i^Xi)/ ^i^x^, where is the shortest wavelength and Xf is the longest wavelength in 
the band. 

the disk between 7 and 32 i?*. With this range of Rq and 
assuming AR ^ 0.1 Rq, the corresponding fiUing factors 
for the 1400 K component range from ~ 10 to 300. 
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Figure 13. Ratios of excess flux on the blue side of the Paschen 
jump (0.815 /im) to that on the red side of the Paschen jump 
(0.915 fim) versus Ey- 

projected area of the star, which can be less than one 
(smaller than the star) or greater than one (larger than 
the star). Typical filling factors for hot accretion spots 
are less than 1% in most cases b ut can be up to 5% in a 
few o f the most active accretors ()CG98l : iGullbring et al.l 
[2000). If a zone in the disk is treated as a face-on, flat 
annulus, its filling factor relative to the star has a simple 
expression. In this case, / = 2RoAR, where Rq is the 
distance from the origin to the center of the annulus, and 
AR is the width of the annulus, both in units of the stel- 
lar radius. For example, an annulus of width 1 i?* at a 
distance of 10 R^, has / = 20, i.e., its are a is 20 time s the 
projected area of a star with radius 1 i?* . IMCHD03I mod- 
eled SpeX 2-5 /im spectra of nine CTTS and found that 
their excess emission over this wavelength range was well 
described by T ~ 1400 K dust from sublimation radii in 



5.1. Constraints from Ey , Ey , an 



d E 



K 



To illustrate the inadequacy of only two regions with 
Thot = 8000 K and Tcooi = 1400 K to account for the ob- 
served lYJ excess, we first look at the median excesses at 
the extremes and center of our wavelength coverage, Ey , 
Ey, and Ek (Fig. [TT|) . before turning to the full shape 
of the excess. In the upper row of Figure [T5l we compare 
the observed excess ratios at pairs of wavelengths to those 
from a two-component model with blackbody tempera- 
tures Thot = 8000 K and Tcooi = 1400 K. The observed 
median relation for Ey /Ek is well described by a com- 
bination of these two temperature regimes if the filling 
factor of the cool component is 1240 times that of the 
hot component. However, the required filling factors for 
the hot region fall between 0.015 (LkCa 8) and 0.22 (DR 
Tau), in contrast to typical va lues of 0.001-0.05 inferred 
from accretion shock models (|CG98[ ). Furthermore, the 
observed ratios Ey j Ey and Ey j Ek are on average twice 
as large as the best two-component model predicts. 

The unexpectedly large Ey measurements suggest that 
a third component, with a temperature intermediate to 
an accretion spot with Thot = 8000 K and sublimating 
dust with Tcooi = 1400 K, must contribute to the ob- 
served excess. We can roughly constrain the allowed 
range of temperatures Tint for this region by compar- 
ing the Ey , Ey, and Ek relations to flux ratios from 
a single-temperature thermal radiator as shown in the 
lower row of Figure [TS] We see that Ey / Ey is too large 
for 8000 K, corresponding instead to a thermal radia- 
tor of 5000 K, and Ey/EK is far too large for 1400 K, 
corresponding instead to a thermal radiator of 2200 K. 
We infer from this that the temperature(s) of the third 
component must lie between 2200 and 5000 K. 

We explore whether a model with three blackbody 
components, characterized by three temperatures and 
three filling factors, can explain the observed excesses. 
In order to constrain parameter space, we specify two 
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Figure 14. Comparison of veiling-corrected equivalent widths of three strong emission lines to the median excesses Ey, Ey , and Ek 
(Table[6ll. The correlation coefficients appear in the upper left of ea ch p anel. All line/continuum pairs are simultaneous, either from SpeX 
or echelle data, except for Br7 (SpcX) and Ey (echelle). As in Fig. 1111 Ey for DL Tau is adjusted to account for variability. 



temperatures from the outset, Thot = 8000 K and Tcooi = 
1400 K, and search for combmations of Ti„t, /hot, /int, 
and /cool that fulfill the observed Ey / Ey and Ey/Ex 
relations. If, for example, we specify Tint, then unique 
ratios of filling factors are required among the three tem- 
perature regions. This is true because the ratio of excess 
fluxes at two wavelengths can be written as 



(Thot) + ^Ai (Tint) /int //hot + -Bai (Tcool) /cool/ /hot 
Bx^ (Thot) + Bx^ (Tint) /int //hot + -^Aa (Tcool) /cool/ /hot ' 

where Bxi{Tj) is the Planck function for temperature 
Tj evaluated at wavelength A^. A system of two such 
equations (e.g., for Ey/Ey and Ey/Ex), can then be 
solved to find expressions for /int //hot and /int //cool that 
depend only on the three assumed temperatures and the 
observed median excess ratios. 



We choose two values for Tnt as an illustration. In 
Case A, Ti„t = 5000 K, as might correspond to lower- 
temperature annuli surrounding the hot accretion spots 
on the stellar surface. In Case B, Tint — 2500 K, as 
might correspond to disk gas inside the dust sublimation 
radius. Given the median observed ratios Ey / Ey = 
0.50 and Ey/Ex = 1.45, the ratios of filling factors 
will be /hot//int//cooi = 1/30/6100 for Case A and 
/hot//i„t//cooi = 1/62/910 for Case B. Although the fill- 
ing factors in these models remain in fixed ratio to one 
another, the magnitudes of the observed excesses require 
that the magnitudes of the filling factors span a range 
of about a factor of 25. The resulting values if the pho- 
tospheric temperature is 4000 K are illustrated in Fig- 
ure [16] and are listed in Table [T] The implied sizes for 
the hot and cool components are reasonable for the phys- 
ical regimes attributed to them in the above scenarios. 
In both cases, the hot-spot filling factors are mostly ^ 1 
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Figure 15. Comparison of observed median excesses E-y , Ey, and Ek to those of blackbodies. Top row: Excesses from two components 
with temperatures 8000 K and 1400 K and a ratio of filling factors of 1240. While this two-component fit can describe Ev/Ek, the 
corresponding ratios Ey / Ey and Ey / Ejf are smaller than observed by more than a factor of two. Bottom row: Correspondence of 
single-temperature blackbodies to excess ratios at different pairs of wavelengths. All excesses are in units of the star's photospheric flux at 
0.8 ^lra. 

with a maximum of 0.23, and the cool dust filling factors 
range from 8 to 280, in line with the crude assumption of 
a dust sublimation region at a distance i?o between 7 and 
32 i?, from the star and with width 0.1 i?o- The cor- ^ Cool 
responding size range for the intermediate-temperature 
component in Case A (Tint — 5000 K) is /int = 0.05 
to 1.4 and in Case B (Tint = 2500K) is /i„t = 0.6 to 
14. For the low-excess stars, if Tint ~ 5000 K, then /int 
would be small enough to correspond to warm annuli 
around accretion hot spots on the stellar surface, but for 

the high-excess stars, this is an unphysical interpretation o.OOi 0.0 1 

since /int can be > 1. In contrast, if Tnt ~ 2500 K, then 
/int is always large in comparison to the star, but smaller 
than a ring at the dust sublimation radius. 

While a single blackbody of Tint with size in exact pro- 
portion to that of accretion spots on the photosphere is 
likely an oversimplified description of the origin of the 
lY J excess for all the CTTS in our sample, we infer 
that all three components have sizes that scale with the 
magnitude of the excess. There is already evidence for 
this in models of accretion shocks, where filling factors 
up to 5% are invoked to explain the largest blue contin- 
uum excesses ( Gullbring ct al. 2000), and in the derived 
dust sublimation radii, w hich are largest for CTTS with 
high disk accretion rates (jEisner et al.ll2010l ). It is sur- 
prising that, for temperatures between 2500 and 5000 K, 
the requisite size of the intermediate-temperature com- 
ponent is comparable to or exceeds the stellar surface 
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Figure 16. Filling factors, in units of the projected area of 
the star, required to describe the observed excesses in a three- 
component model. The solid (Case A; Tjut = 5000 K) and dashed 
(Case B; Tint = 2500 K) lines show the ranges for all stars (Ta- 
ble [Tj. Diamonds mark the individual fill ing factors for the four 
stars with SEDs modeled in § 5.2 (Fig. [TSll. 

area for those stars with large lYJ excesses. In the next 
subsection we go beyond the simple median ratios of Ey, 
Ey, and Ek, and we experiment with fitting the full SED 
of the excess emission. 

5.2. SED Fitting 

Given the need for at least three temperature compo- 
nents with sizes that increase with the magnitude of the 
excess to account for the median ratios among Ey, Ey, 
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Table 7 

Filling Factors Yielding Median Excess Ratios for Fixed Ti^t 



Case /int//hot /cool//hot /cool//int Range of /hot Range of /int Range of /cool 

A: Tint = 5000 K 30 6100 200 0.002-0.05 0.05-1.4 11-280 

B: Tint = 2500 K 62 910 15 0.009-0.23 0.56-14 8.2-210 



and Ek, we attempt here to reproduce the full shape of 
the excess emission XEx. We continue to use the Case A 
and Case B scenarios described in the previous section, 
with Thot = 8000 K, Tcooi = 1400 K, and Tint = 5000 K 
(Case A) or 2500 K (Case B). 

A comparison of calculated spectra from Case A and 
Case B to the median excesses at seven wavelength in- 
tervals approximating the R, I, Y, J, H, and K 
bandpasses (Table [6] and Fig. [TO)) is shown in Figure [iTl 
Either scenario accounts reasonably well for the general 
behavior of the observed excesses, but Tint — 5000 K is 
a better match to the excesses between V and J, and 
Tint — 2500 K is a better match between J and K. 

Further comparison of Cases A and B to the observed 
excesses can be made by dropping the requirement of 
a uniform ratio of filling factors for all stars and let- 
ting the filling factors for each region match the rela- 
tions among Ey, Ey , and Ek for each individual star, 
which can differ by factors of several from the median 
ratios. We show this in Figure [18] for four representa- 
tive sources: two high-excess stars (DR Tan and CW 
Tau) and two moderate-excess stars (DK Tau A and BP 
Tau). We compare their full observed \E\ spectra de- 
rived from the echelle and SpeX data to both Case A 
and Case B models. The filling factors chosen for each 
star are listed in Table |5] and marked in Figure [Tni In 
addition to the summed fit, we also show the individual 
contribution from each component. Again we see that, 
although these simple three-component models give rea- 
sonable fits to XE\ , the warmer Case A better represents 
the shorter- wavelength excess, while the cooler Case B 
better represents the J through K bands. We conclude 
that, although Cases A and B both describe the observed 
excesses reasonably well, the source of the excess emis- 
sion in this wavelength regime is more likely to arise from 
regions with a range of temperatures between 2200 and 
5000 K. 

5.3. Residual Spectrum and Luminosity 
of the lYJ Excess 

To uniquely determine the excess emission spectrum 
that arises solely from the region of intermediate temper- 
ature, one must specify the contributions from the hot 
and cool components, which would require even greater 
wavelength coverage than we present here. We can, how- 
ever, constrain the spectral shape and luminosity of the 
intermediate-temperature component by estimating the 
contributions from the hot and cool components and sub- 
tracting these contributions to generate a residual excess. 
We show such residual spectra in the first two columns 
of Figure [19] for the four stars whose excesses were fit 
with Case A and B scenarios in Figure [18] For clar- 
ity, only SpeX data were used where they overlap with 
echelle data, and the XEx curves were median-filtered 
and rebinned to a sparse wavelength grid in order to re- 
move emission lines, placing the focus on the overall con- 



tinuum shape. For each star, the residual (solid line) is 
compared to the blackbody (dotted line) for the assumed 
Tint of Case A or B. We find good agreement between the 
residual spectrum and the corresponding blackbody. 

An alternate approach is to define a minimum resid- 
ual emission by assigning all of the excess at 0.48 /im to 
an accretion spot with a blackbody spectrum of Thot = 
8000 K, assigning all of the excess at 2.4 /xm to dust with 
a blackbody spectrum of Tcooi = 1400 K, and then re- 
moving their contributions to the observed XE\ between 
0.48 and 2.4 /im. By comparison, in Case A the flux at 
0.48 fim has a substantial contribution from both the in- 
termediate and hot components, and in Case B the flux 
at 2.4 ^m has a substantial contribution from both the 
intermediate and cool components (see Fig. [TS]). These 
minimum residuals are shown in the last column of Fig- 
ure]!!] (solid lines) and are compared to blackbodies (dot- 
ted lines). The blackbody spectra are in good agreement 
with the minimum residuals, with temperatures interme- 
diate to Cases A and B: 3900 K for DR Tau, 3470 K for 
CW Tau, 3600 K for DK Tau, and 3980 K for BP Tau. 

The luminosities corresponding to these three ap- 
proaches for evaluating the residual emission from the 
intermediate-temperature component are shown in each 
panel of Figure 1191 Luminosities are calculated by nu- 
merical integration of the excess emission spectra, ini- 
tially expressed in units of L* as ratios of excess luminosi- 
ties to photospheric luminosities. To convert to Lq, these 
are multiplied by literature estimates of the stellar lumi- 
nosities (TableU]). For each of the four stars, the range in 
luminosities from the three approaches is less than a fac- 
tor of two. It is obviously of interest to compare the esti- 
mated luminosities of the intermediate-temperature com- 
ponents to previously published accretion luminosities. 
For all 14 CTTS in our sample that have both echelle 
and SpeX spectra, we do this by estimating both the 
minimum and maximum luminosity of the intermediate 
component from the observed excess. For the minimum 
luminosity we use the approach just described for all 
stars, forcing the intermediate-temperature contribution 
at the shortest and longest wavelengths to zero. For the 
maximum luminosity we attribute the entire excess be- 
tween 0.48 and 2.4 /zm to the intermediate-temperature 
component. Again, the luminosities are first found in 
photospheric units and then converted to physical units 
using published values of the photospheric luminosities. 
The minimum and maximum estimates, along with the 
stellar and accretion luminosities from the literature, are 
listed for each star in Table 15] and compared graphically 
in Figure [20] The minimum luminosities run from 0.01 
to 2.4 Lq (0.03 to 1.2 L*), while the maximum luminosi- 
ties run from 0.09 to 7.8 Lq (0.22 to 4.2 L*). Literature 
accretion an d phot ospheric luminosities were taken from 
IGHBC98I or ICG98I whenever possible (aU but two stars) 
in order to ensure a uniform approach. Figure shows 
that the minimum and maximum luminosity estimates of 
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Figure 17. Comparison of -Baa for Case A (Tint = 5000 K) and Case B (Tjnt = 2500 K) to observations (gray lines, see Fig. 1101 for 
star identification). The strengths of the model excesses are characterized by the corresponding veiling at the shortest wavelength of our 
dataset, = 0.3 to 6.3. Both three-component fits are reasonable approximations to the observed spectra, although Case A is a better fit 
at short wavelengths and Case B is a better fit at long wavelengths. 



Table 8 

Three-Component Fit Parameters for Selected CTTS 



Object 


Case 


/int //hot 


/cool / /hot 


/cool/ /int 


/hot 


/int 


/cool 


DR Tau . . 


A 


12.5 


3310 


265 


0.080 


1.00 


265 


CW Tau . 


A 


23.0 


10000 


435 


0.027 


0.62 


270 


DK Tau A 


A 


18.8 


6130 


327 


0.008 


0.15 


49 


BP Tau . . 


A 


17.5 


3000 


171 


0.008 


0.14 


24 


DR Tau . . 


B 


55.0 


1050 


19.1 


0.20 


11.0 


210 


CW Tau . 


B 


80.9 


1910 


23.6 


0.11 


8.9 


210 


DK Tau A 


B 


70.3 


1370 


19.5 


0.03 


1.9 


37 


BP Tau . . 


B 


60.0 


640 


10.7 


0.03 


1.5 


16 



the lYJ excess bracket the accretion luminosity deter- 
mined from the blue excess, suggesting that the actual 
luminosity of this component is comparable to the pre- 
viously derived accretion luminosity for each star. The 
implication is that total accretion luminosities for most 
CTTS are a factor of two higher than previously esti- 
mated from modeling the blue excess emission. 

In summary, from our very simple modeling of the ex- 
cess emission between 0.48 and 2.4 fim, we find that three 
temperature components whose filling factors increase 
with increasing total excess account reasonably well for 
the observations. The temperature structure of the in- 
termediate component is not well specified, but gas with 
one or more temperatures between 2200 and 5000 K is 
likely. It seems inescapable that for stars with the largest 
lYJ excesses (highest disk accretion rates), the requisite 



filling factor is comparable to the surface area of the star 
if Tint is 5000 K and is up to an order of magnitude higher 
if Tint is only 2500 K. This is unreasonably large to at- 
tribute to warm annuli around spots in the shock-heated 
photosphere and points to an origin in the disk, accretion 
flow, or wind. The luminosity from this region is compa- 
rable to previous estimates of the accretion luminosity. 

6. DISCUSSION 

We have demonstrated the presence of an excess con- 
tinuum in accreting T Tauri stars between 0.48 and 2.4 
fim with an intensity that scales with other indicators 
of disk accretion, such as the blue excess continuum and 
atomic emission lines. The similarity in the gross spectral 
characteristics of the lYJ excess among a diverse sam- 
ple of stars makes it unlikely that contamination from 
random factors, such as large cool spots or undetected 
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Figure 18. Comparison of tliree-componcnt blackbody fits for Case A (Tint = 5000 K; left) and Case B (Tint = 2500 K; right) to tlie 
emission excess \Ex (solid lines, from Fig. |9ll for four stars. Dotted lines show the three components of each model, and dashed lines show 
their sum. The required filling factors for each component are labeled (also see Table [S]!. 



companions, contributes much to the derived excess. We 
discuss some imphcations of the lYJ excess, which has 
a temperature between that attributed to ^ 8000 K hot 
spots in the shock-heated photosphere (the dominant 
contributor to the excess shortward of 0.5 /im; CG98) 
and a raised rim of ~ 1400 K dust at the dust subh- 
mation radius in the disk (the dominant contributor to 
the excess from 2 to 4.8 /xm; IMCHD03I '). Our simple 
modehng suggests that the intermediate component has 
a temperature between 2200 K and 5000 K, and it hkely 
has a range of temperatures between these hmits. One 



issue to consider is the effect of this newly recognized 
contribution to CTTS excess continuum emission on the 
determination of disk accretion rates. Another is the 
absence of any spectral region where broadband colors 
sample the young stellar photosphere, which makes ex- 
tinction determinations from broadband colors difficult. 
Most exciting is the potential of identifying a heretofore 
unrecognized aspect of the star-disk interaction in accret- 
ing systems. 
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Figure 19. Residual XEx for the intermediate-temperature component with three different assumptions for four rep rese ntative stars. The 
Case A {left) and Case B (center) residuals are derived by subtracting the hot and cool excess contributions in Fig. 1181 from the observed 
excess. Minimum-luminosity residuals are in the right column, for the extreme assumption that the excess at 0.48 fim is entirely from the 
hot component and at 2.4 fim is entirely from the cool component (see Table[9ll. Dotted lines correspond to blackbodies of Tint = 5000 K 
(left), Tint = 2500 K (center), and best fits (right), where the latter are between 3500 and 4000 K. 



6.1. Effect of the lYJ Excess on Derived 
Disk Accretion Rates 

Although wc have identified a new source of continuum 
emission between 0.48 and 2.4 fim, we cannot uniquely 
determine its properties due to spectral overlap with 
other known sources of excess emission. However, rea- 
sonable upper and lower limits indicate the excess lu- 
minosity in this wavelength range is comparable to pre- 
viously determined accretion luminosities. Provided the 



excess is not simply reprocessing radiation from accretion 
hot spots, this indicates that a significant amount of ac- 
cretion energy is being dissipated in the region we call 
the intermediate-temperature component. It is not im- 
mediately obvious what the implication is for estimates 
of disk accretion rates. When disk accretion rates are 
determined from spectrophotometry blueward of 0.5 /im, 
this em ission is well represented by accretion shock mod- 
els (jCG98l ) where the underlying assumption is that the 
excess luminosity derives from kinetic energy liberated 
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Table 9 

Minimum and Maximum Excess Luminosities 



Object 


v-^w / 


/ J ^r^r- i 1 ] 

-^acc V W / 




r . (T \h 


/ Jrrt ; 11 I 1 Jl^i 1 


r If Y 


T (T r-:\ 


AA Tau . . 


0.71 


0.03 


2,2 


0.08 


0.06 


0.14 


0.10 


AS 353A . 


3.72 




3 


0.66 


2.44 


2.10 


7.82 


BM And . 


7.24 




5 


0.07 


0.50 


0.31 


2.28 


BP Tau . . 


0.93 


0.18 


2,2 


0.16 


0.15 


0.49 


0.46 


CW Tau. . 


1.35 


1.41 


4,1 


0.35 


0.47 


1.34 


1.81 


CY Tau . . 


0.46 


0.04 


2,2 


0.06 


0.03 


0.28 


0.13 


DF Tau . . 


1.97 


0.36 


2,2 


0.24 


0.47 


0.67 


1.31 


DG Tau A 


1.74 


4.96 


3,1 


0.23 


0.39 


0.84 


1.46 


DK Tau A 


1.45 


0.17 


2,2 


0.20 


0.29 


0.68 


0.98 


DL Tau . . 


0.68 


0.61 


3,1 


0.75 


0.71 


1.72 


1.17 


DO Tau . . 


1.01 


0.60 


2,2 


0.22 


0.22 


1.13 


1.14 


DR Tau . . 


1.74 


2.74 


3,1 


1.15 


2.00 


4.17 


7.26 


HN Tau A 


0.19 


0.04 


2,2 


0.12 


0.02 


1.37 


0.26 


LkCa 8 . . . 


0.41 


0.01 


2,2 


0.03 


0.01 


0.22 


0.09 



Re fer ences. — (1) ICG98I : (2) EHBCSI f3') IHartigan et al.|[T99a: C4'l IKenvon fc Hartmaiiiil 
[T995I : (b) IRostoDchina|[T999l . 

References for stellar and accretion luminosities. 
^ Intermediate-component luminosity if the excess emission at 0.48 fim is due entirely to hot 
spots and the excess emission at 2 A fim is due entirely to warm dust. 

Total excess luminosity from 0.48 to 2.4 fim. 

the excess will come from the intermediate-temperature 
component. In order to correctly determine CTTS accre- 
tion rates based on line veilings, the total accretion lumi- 
nosity from all temperature components would first need 
to be found from absolutely simultaneous spectropho- 
tometric data extending from shortward of the Balmcr 
jump in the U band to longward of the gas continuum 
excess, probably out to the K band. Evaluation of red- 
dening with a good grid of both subdwarf and dwarf stan- 
dards through this full wavelength region would be nec- 
essary to understand the tendency to find larger Ay at 
longer wavelengths in both WTTS and CTTS (see next 
section) and thus to appropriately correct the observed 
emission for extinction. Then to convert line- veiling mea- 
surements to accretion luminosities, the relationship be- 
tween veiling observed over a small wavelength interval 
and the total accretion luminosity must be established 
in order to know what kind of bolometric correction to 
apply. Finally, a model that includes all sources of excess 
luminosity would need to be constructed so that accre- 
tion luminosities can be turned into disk accretion rates. 

At present, we can conclude only that the excess lYJ 
emission has a luminosity that is comparable to previ- 
ously derived accretion luminosities and that accretion 
rates based on line veiling longward of 0.5 /xm are sus- 
pect. An additional problem is that, for many CTTS, 
photospheric properties such as luminosity, radius, and 
mass will need to be re-evaluated in light of the fact that 
much of the flux between 0.5 and 1 /im previously at- 
tributed to the photosphere is actually excess emission. 
Flux-calibrated spectrophotometry from the UV to the 
near IR in combination with sophisticated models of all 
contributors to the excess emission will be needed to im- 
prove our understanding of accretion luminosities and 
disk accretion rates. 

6.2. Effect of the lYJ Excess on Extinction Estimates 

The demonstration that CTTS have excess contin- 
uum emission in the region intermediate to the well 
studied blue and infrared excesses, including the /, Y, 
and J bands, makes it clear that there are no wave- 
length regions where photometric colors sample the bare 
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Figure 20. Comparison of minimum and maximum estimates for 
the luminosity of the intermediate-temperature component giving 
rise to the lY J excess for each star to previously published ac- 
cretion luminosities. The maximum assumes that all of the excess 
between 0.48 and 2.4 fiva comes from the intermediate-temperature 
component, and the minimum assumes that all of the excess at 0.48 
^m comes from an 8000 K accretion shock, while all of the excess 
at 2.4 /jm comes from 1400 K warm dust (also see Table|9ll. 

by material free-falling to the stellar surface in accre- 
tion columns. Thus accretion rates determined under 
these assumptions should be reasonable estimates of the 
amount of matter actually landing on the stellar surface 
if the TY J excess arises from a region other than the ac- 
cretion shock. However, if some of the lY J excess arises 
from annuli of warm gas surrounding a central, small ac- 
cretion hot spot at the magnetic footpoint, then accretion 
rates will need to be revised upward accordingly. 

More troubling is the commonly used technique of esti- 
mating disk accretion rates from large samples of CTTS 
based on line- veiling measurements in the Paschen con- 
tinuum longward of 0.5 /zm. Since this approach relies on 
a bolometric correction appropriate for a hot accretion 
shock, it will overestimate accretion luminosities, pro- 
gressively more so at longer wavelengths, since much of 
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Figure 21. The difference between the apparent reddening Ay^'^'^ and the true reddening yl^''^"'^ as a function of color if observed CTTS 
colors are assumed to be photospheric. The underlying photosphere is 4000 K and the excess continuum is the sum of three components 
corresponding to Cases A (black) and B (re d). The magnitude of the excess is identified by the veiling at V, with ry from 0.3 to 6.3, 
corresponding to the SEDs shown in Fig. 1171 



photosphere. Thus, using observed colors to determine 
Ay is subject to error. Many compilations of extinc- 
tions for CTTS use red colors t o deriv e Ay, such as 
V-R (IKenvon fc Hartmandll995l) . V-I (IWhite fc Ghez 
200 l . -R- J (|Eiroa et al.ll2002D . or 7 - J (|McClure et al. 
2010[) . To assess the effect of the lYJ excess on deter- 
mining Ay from observed colors, we synthesized SEDs 
from 0.5 to 2 /im that are composed of a 4000 K stellar 
photosphere with varying amounts of continuum excess 
and subjected to intervening extinction. The continuum 
excess consists of three temperature components corre- 
sponding to the Case A and Case B of the previous sec- 



tion, with fixed Thot and Tc 



and Ti, 



2500 or 5000 



K. The magnitude of the excess is traced by the veiling at 
with ry ranging from 0.3 to 6.3. The resulting differ- 
ence between the apparent and true Ay as a function of 
color is shown in Figure [H] A failure to correctly account 
for the lYJ excess emission leads to an underestimate of 
Ay from all colors except I — J, which are often not ob- 
served simultaneously. Although there is excess emission 
at / and J, the intermediate components in Cases A and 
B have temperatures similar to the star, and thus the 
I — J color, which primarily samples the intermediate 
component, happens to be nearly photospheric. On the 
other hand, the underestimate of Ay for heavily veiled 
stars can be as large as 2 magnitudes when using V — R 
and 1 magnitude when using R — I. (The difference be- 
tween the true and derived Ay is equal to the difference 
between the V magnitude of the bare photosphere and 
the V magnitude of the excess-plus-photosphere system; 
thus, it depends only on the properties of the excess and 



is independent of the true Ay.) 

An additional issue regarding extinction measurements 
in CTTS is the puzzling tendency to find Ay from 0.8 
to 2.4 /im spectra that consistently exceed previous es- 
timates. This could result from a number of effects that 
may all be operating in concert, including incorrect as- 
sessment of Ay from red colors due to a large lYJ ex- 
cess, differences in surface gravity or temperature be- 
tween the CTTS and dwarf spectral templates, and dis- 
tortion of the spectrum by large cool spots or undetected 
close companions. (We consider the last two options un- 
likely because of the similarity in the gross character- 
istics of the lYJ excess in our 16 CTTS.) This effect 
may be related to the unexplained behavior in the SEDs 
of some WTTS, which have increasingly anomalous col- 
ors for their spectral ty pes toward longer wavelengths 
(IGHBC98I: ICullbring et al. 19981^. We tested the SpeX 
spectrum of the WTTS V819 Tau to see if this behavior 
could be explained by an lYJ excess that arises from a 
blackbody source with T = 2500, 3500, or 5000 K and is 
weak enough to produce negligible veiling shortward of 
0.6 /im. We find that such a weak excess does generate a 
higher Ay with increasing wavelength, but only by a few 
tenths of a magnitude, much smaller than observed. We 
reiterate that whatever gives rise to the large Ay deter- 
mined from the SpeX spectra does not significantly affect 
our characterization of the lYJ excess, as described in 
§ 4, but it does indicate that extinctions to CTTS and 
WTTS need to be revisited in a more thorough way than 
can be done here. 
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6.3. A New Component in Accretion Disk Systems? 

Identifying the source of the lYJ excess continuum re- 
quires that we define its SED in order to determine its 
temperature, density, opacity, and filling factor. Since at 
the shortest wavelengths in our data there is a contribu- 
tion to the excess from the hot accretion spots, and at 
the longest wavelengths there is a contribution from the 
warm dust rim in the disk, we cannot extract the residual 
SED arising solely from the additional red emission. It 
does appear to be broad and relatively smooth, as would 
be expected from blackbody radiation with a tempera- 
ture (or temperatures) between 2200 and 5000 K. Pos- 
sible sources for the red continuum emission are warm 
annuli around accretion hotspots, gas in the disk in- 
side the dust sublimation radius, and gas in the accre- 
tion flow or wind. The presence of multiple temperature 
components in the shock-heated photosphere would not 
be surprising; in fact, this is predicted by 3D magneto- 
hydrodynamic numerical simulations of accretion in the 
presence of dipolar arid more cornple x magnetic fields 
(jRomanova et al.l 120041 : iLong et all 12008) . This model 
may apply to stars with small lYJ excesses, provided 
that about half of the accretion luminosity is radiated 
from these cooler regions. However, for stars with large 
lYJ excesses, the required filling factors exceed the size 
of the star, making this an improbable scenario. The 
areas of the accretion flow and wind are larger than the 
surface area of the star, making these plausible locations 
for the lYJ excess, although the requisite temperature 
is much cooler than the 10,000 K that is invoked for 
accretion streams as a source of hydrogen line emission 
(jMuzerolle et al.|[200l . 

The possibility that the lYJ excess is due to dust-free 
gas in the disk inward of the raised rim at the dust sub- 
limation radius deserves consideration. Theoretical con- 
siderations suggest that the dust-free inner gaseous disk 
should be optically thin with strong molecular emission 
(jMuzerolle et al. 2004), but this conclusion is called into 
question by IDullemond fc Monnierl (2010). These au- 
thors note that there appears to be a deficit of molecules 
in the dust-free inner disk, which only rarely shows CO 
overtone or H2O emission (jNaiita et al.ll200"7l) . and they 
suggest that molecules are destroyed in the inner gaseous 
disk, possibly from proximity to photodissociating radi- 
ation. They also present evidence that the inner gaseous 
disk may be optically thicker in CTTS than in interme- 
diate mass Herbig Ae/Be stars, such as the behavior of 
linear polarization across the Ha line, which appears to 
be sc attered off a rotating inner gaseous disk (iVink et aLl 
[2005f ). 

Additional evidence for disk gas inside the dust 
sublimation radius comes from near-infrared spectro- 
interferometry in the K band. Both spectroscopy and 
near-infrared interferometry indicate that the dust subli- 
mation zone is at a rad i al distance in the di s k of ~ 0.1 AU 
for C TTS (MCHD03" 'Eisner et al.' 2001 lAkeson et afl 
120051 iMillan-Gabet e t al. 2007). Near- infrared spectro- 
interferomet ry at a resolution of R = 230 across the 
K band by lEisner et al.l 120091 (and lEisner et al] 120071 
& iTannirkulam et al.ll2008al for higher-mass AeBe stars) 
shows that the angular size of the near-infrared emission 
decreases with decreasing wavelength, indicating that the 
ring of warm dust is accompanied by warm, presumably 



gaseous material within the dust sublimation radius. For 
higher-mass stars than we have investigated here, the 
possibility that this excess may arise from non-blac kbody 
small d ust grains has also been considered (Benist v et al. 
20T0a"P). Detailed modeling (e.g., bv ITannirkulam et al. 
2008b) indicates that for gas opacity primarily due to 
free-free H and H~ emission, the gas is partially opti- 
cally thin with continuum emission primarily from ma- 
terial hotter than 2000-3000 K but cooler than 5000 K. 
It is possible that the excess continuum in our broader 
spectral coverage comes at least in part from this gaseous 
region in the inner disk. 

The inner gaseous disk, lying between the dust subli- 
mation zone and the truncation of the disk by the stel- 
lar magnetosphere, is not yet well characterized but is a 
likely place for accretion energy to be dissipated. The 
gas that accretes onto the star must pass through this 
region, and magneto-centrifugual disk winds could be 
launched here. If the lYJ continuum emission arises 
at least in part in this region, line emission would also 
be expected. Alth ough the region appears to be de- 
pleted of molecules (jPuUemond fc Monnieril201 0^). some 
of the abundant atomic emission lines in CTTS spec- 
tra may arise here. Near-infrared spectro-interferometry 
shows that the Brackett emission in CTTS is unre- 
solved ()Eisner et al.l I2010D . requiring formation inside 
the resolved dust sublimation region, and in this study 
we found a tight correlation between the lYJ excess 
and the Paschen and Brackett line fiuxes. Although 
the kinematic structure of these lines points to an ori- 
gin in the acc retion flow and inner wind rather than 
the inner disk (iMuzerolle e t al.l l 2001l: iFolha fc EmersonI 
2OOIHK wan fc F ischer 2010), the unexpected finding of 
Edwards et al.l (|2006i) that the structure of the Pa7 line 
profile depends on the 1 /im veiling, with the broadest 
lines in stars with the largest veiling, makes it clear that 
the origin of these lines is not yet fully understood. Thus, 
a contribution from the inner disk cannot be ruled out. 
More broadly, the complexity of the rich emission line 
spectrum of CTTS in the optical and near infrared is 
not fully explored and may yet reveal li nes that are pri- 
maril y formed in the inner gaseous disk (jBeristain et al.l 

mil). 

Ongoing efforts to definitively identify the physical 
properties of the inner region of accretion disks should 
yield an enriched understanding of accretion in CTTS. 
The lYJ excess is of particular interest in the CTTS that 
are the most heavily veiled, such as DR Tan, DC Tau A, 
RW Aur A, and AS 353A. Such CTTS are borderline 
class I/II young stellar objects, with large infrared ex- 
cesses, large accretion rates, and spatially resolved micro- 
jets, and their large lYJ excesses can provide a view of 
the star-disk interaction region in a phase of protostellar 
evolution that is usually heavily obscured in the optical 
and near infrared. 

7. CONCLUSIONS 

We have characterized the excess continuum emission 
in accreting T Tauri stars between 0.48 and 2.4 /im. At 
1 /Ltm the magnitude of the excess ranges from 0.1 to 3 
times the photospheric flux, and it scales with the excess 
at shorter and longer wavelengths. We interpret this ex- 
cess as arising from a source of continuum emission in 
addition to small hot (T ~ 8000 K) accretion spots on 
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the stellar surface at the base of accretion columns and 
large rings of warm (T ^ 1400 K) dust at the sublima- 
tion radius in the disk. The broad, smooth distribution 
of the lY J excess suggests it comes from a third compo- 
nent with a temperature between 2200 K and 5000 K and 
roughly blackbody emission characteristics. We cannot 
pin down a unique combination of temperature and fill- 
ing factor to account for the emission, but among stars 
with high lY J excess, the filling factor is comparable to 
the surface area of the star if the temperature is 5000 K 
or more than an order of magnitude larger if the tem- 
perature is 2500 K. Possible sources of the emission in- 
clude warm annuli surrounding hot accretion spots in the 
shock-heated photosphere, accreting gas in funnel flows, 
or disk gas inside the dust sublimation radius. The lu- 
minosity from this region is comparable to the accretion 
luminosity found from ultraviolet and blue excesses, sug- 
gesting that accretion luminosities have been underesti- 
mated by a factor of about two. The implication for disk 
accretion rates is more complicated, but a comprehen- 
sive campaign to obtain simultaneous spectrophotometry 
from the ultraviolet to the near infrared for a large sam- 
ple of CTTS would provide further insight. Whatever its 
source, this excess emission contains a large fraction of 
the accretion energy released in the vicinity of the star. 
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